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CROSS-REFERENCE TO RELATED APPLICATIONS 

[001] This application claims priority of Korean Application No. 

10-2003-0056189, filed on August 13, 2003, the disclosure of which is 
incorporated fully herein by reference. 



FIELD OF THE INVENTION 

[002] Generally, the present invention relates to a vehicle. More particularly, 

the present invention relates to a suspension of a vehicle and a method for controlling 
the same. 

BACKGROUND OF THE INVENTION 

[003] A suspension of a vehicle responds to various road impacts encountered 

while the vehicle is being driven, and controls dynamics of the vehicle under 
acceleration/deceleration and turns. One important role of such a suspension is to reduce 
transmission of noise and vibration into a vehicle body. 

[004] In a typical suspension of a vehicle, a spring and a damper (also called a 

shock absorber) is mounted to each wheel. A spring coefficient (also called a spring rate 
or a spring constant) of the spring, a damping coefficient of the damper, and mounting 
positions of the springs and dampers are determined in the designing process of the 
suspension. Therefore, a great deal of research has been undertaken to enhance 
performance of a suspension by efficiently determining such design factors. 
[005] However, the amount of possible enhancement of the performance of a 

suspension by simply changing such design factors is limited. So recently, new research 
has been undertaken to actively eliminate vibration of the suspension by utilizing a 
plurality of excitation actuators mounted to a suspension system. Such excitation 
actuators may be disposed with each spring/damper for each wheel. 
[006] For example, through a process called direct velocity feedback control, 

vibration energy may be dissipated by increasing the damping effect in a vehicle 
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suspension system that has multiple degrees of freedom. Such a technique shows a merit 
in that asymptotic stability is ensured. 

[007] However, in order to obtain asymptotic stability by such direct velocity 

feedback control, many sensors and excitation actuators are required. This is 
problematic because equations of motion of a suspension become inter-coupled when 
modeling the vehicle suspension by equations of motion. 

[008] Therefore, recently, research has been undertaken to design a suspension 

such that the equations of motion thereof become non-coupled, or to determine what 
actuating forces should be applied to such a designed suspension. 

[009] The information disclosed in this Background of the Invention section is 

only for enhancement of understanding of the background of the invention and should 
not be taken as an acknowledgement or any form of suggestion that this information 
forms the prior art that is already known to a person skilled in the art. 

SUMMARY OF THE INVENTION 

[0010] An exemplary method for designing a vehicle suspension system 

according to a embodiment of the present invention includes: formalizing the vehicle 
suspension system in accordance with a preferred equation of motion. Such an equation 
is preferably a linear matrix equation having n degrees of freedom (where n is greater 
than 1) and including a damping matrix for a viscous damping. Eigenvectors of a 
stiffness matrix K for the equation are calculated. The eigenvectors are normalized 
with respect to a mass matrix M. A similarity transform matrix S is calculated 
consisting of the normalized eigenvectors. The equation of motion is then normalized 
using the similarity transform matrix *S*. According to such a design process, a 
suspension system designed as such has its equation of motion adaptable for mode 
normalization. 

[0011] In a preferred embodiment of the invention the preferred equation of 

motion is: 

(1) 

Afx(t) + C(x(0 - u(t)) + K(x(t) - 1/(0) = Pf(0 
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[0012] Here, n and p respectively denote the degrees of freedom of the 

suspension system and the number of independent actuators. 

[0013] M, C, and K respectively denote a mass matrix, a damping matrix, 

and a stiffness matrix, each of which is symmetrically nxn. The mass matrix M is a 
positive definite matrix, the damping matrix C is a positive semi-definite matrix, and 
the stiffness matrix K is a positive definite matrix. 

[0014] P denotes an nxp real matrix corresponding to positions of the actuators. 

[0015] -X-CO and ^(0 respectively denote nxl state and disturbance vectors, 

and /(0 denotes a pxl external force (i.e., an actuating force) vector. 

[0016] In a further embodiment, the normalizing of equation 1 establishes a 

k . = a x c . k 

proportional relationship J J between each pair of a spring coefficient J 

Q 

of a j-th spring and a damping coefficient J of a j-th damper corresponding to the j-th 
spring. 

[0017] An exemplary vehicle suspension system according to an embodiment of 

the present invention includes: a plurality of springs; a plurality of dampers each 
corresponding to one of the springs; and a plurality p of actuators for generating control 
force to the suspension system. The suspension system is formalized by the above 
equation 1, and equation 1 is decoupled into n modal equations, wherein equation 1 is a 
linear matrix equation having a number, n, of degrees of freedom, and the linear matrix 
equation includes a damping matrix for viscous damping. 
[0018] In the same way as described above, in a further embodiment, a 

k . = OC XC . 

proportional relationship J J is satisfied between each pair of a spring 

k C 
coefficient j of a j-th spring and a damping coefficient J of a j-th damper 

corresponding to the j-th spring. 

[0019] When the number n and the number p are equal, it is preferable that the 

suspension system further includes: a detecting unit for detecting at least one of the state 

vector and its velocity ; and a controller for controlling the actuators on the 
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basis of the detected at least one of the state vector or its velocity Here, 
the controller controls the actuators by an actuating force of , wherein: 

[0020] Q = S T P fj = ~C S i4i t and x(t) = S %(t) ^ satisfied 

c 

[0021] Si i s a damping coefficient of a sky-hook damper connected to an i-th 

mode. S is a matrix consisting of eigenvectors of the stiffness matrix K and is 
normalized with respect to the mass matrix M. 

[0022] When the number p is less than the number n, it is preferable the 

suspension system further includes: a detecting unit for detecting at least one of the state 

vector and its velocity ; and a controller for controlling the actuators on the 

basis of the detected one of the state vector or its velocity 

x(0 

, wherein the 

controller controls the actuators by an actuating force of that satisfies 

l = -F sl sign(l) = %Qf 

M . Here, Q = S P and x(t) ="S#(0 are 

satisfied. Si is a frictional force of a sky-hook coulomb friction damper connected to 



an i-th mode. S is a matrix consisting of eigenvectors of the stiffness matrix K and 
normalized with respect to the mass matrix M. 

[0023] In more detail, the actuating force preferably satisfies 

VQ^^)^oa^/gi^)^0&..^(gi«J^C|; fj =-F A 
[irfl^gwg)<0^^)<0&..^ignK,)<Q fj=-F B 

with respect to 1 ~\- ^ and J ~ \> ->P. 



is 
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[0024] Here, 17 A is a value in a range of zero (0) to P\ B is a value in a 

range of zero (0) to ^N- F k for k = 1, . . . ,(2 - 2) ig a ya]ue between and F N; 

F F 

and P and N respectively denote a positive maximum force and a negative 
maximum force that a j-th actuator can generate. 

[0025] In further detail, the actuating force preferably satisfies 

PQfignii^O fori=\^ fj=-F A 

ifQ^ign{i)<0fori=\..^ fj=-F B 

Otherwise , / . = 0 

j j 

with respect to ' = \-f l and J = \"*iP. 

[0026] It is further preferable that values of r A and P are equal, and that 

fT J7 
values of B and are equal. 

[0027] A vehicle suspension system controlled by an exemplary method of the 

present invention includes a plurality of dampers and a plurality of actuators, is 
formalized by the above-mentioned equation (1), and is transformed to a second 
decoupled equation. 

[0028] An exemplary method for controlling such a vehicle suspension system 

includes: calculating a velocity vector 

x(t) 

of a state vector of the above 

mentioned equation (1); calculating an actuating force such that the actuating 

force / (0 satisfies /(0 = & P)'' (-C SI )(S T KS)~^ (.S T K)x(.t), the C Si 
being a damping coefficient of a sky-hook damper connected to an i-th mode; and 
actuating the actuators by the calculated actuating force AO. 



5 



[0029] In a further preferred embodiment of the invention, the second decoupled 

equation is: 

(2) 

/|(o + dflgtert l#(o - m)+ a k («o - 7(0) = ho . 

[0030] Here, I is an nxn unit matrix. S is a matrix consisting of eigenvectors 

of the stiffness matrix K and is normalized with respect to the mass matrix M. 

[0031] Q = S T F, f = Of it), x(t) = S %(t) u(t) = Srj(t) 
S T KS = diagico, 2 ] = A K , and S r CS = C = diag^.co, ] are satisfied 

by the matrix S. 

[0032] An exemplary method for controlling such a suspension system may 

include: calculating a velocity vector of a state vector of an equation (1); 

calculating an actuating force fit) such that the actuating force satisfies 
^Q/^i)^0&^ /, =-F A 

ifQjSigK&>0&Q 2J sigK£ 2 )>0&- • Q^igKt) <0, /, =-F x 



ifQ J sign($)<0&Q 2J sigrj(i 2 )<0&- ■ Q nj sign&)>0, /, =-F^_ 2) 

with respect to ' r = \- • and J \*"*P\ and 



actuating the actuators by the calculated actuating force 



fit). 



[0033] Here, r A is a value in a range of zero (0) to P\ r B is a value in a range 

of zero (0) to ^V; for k = \,...^2 — 2) ig a yalue between and and 

F F 

P and respectively denote a positive maximum force and a negative maximum 
force that a j-th actuator can generate. 

[0034] It is preferable that the actuating force satisfies 
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ifQijSignii^O fori=\.../i> f y =-F A 
ifQ,sign{l)<0 fori=\.^ f y =-F B 
Otherwise , f . = 0 

with respect to ' ~\>-fl and J ~"lv?P. 

[0035] It is further preferable that values of r A and P are equal, and values of 

TT 77 
B and JV are equal. 

BRIEF DESCRIPTION OF THE DRAWINGS 

[0036] The accompanying drawings, which are incorporated in and constitute a 

part of the specification, illustrate an embodiment of the invention, and, together with 
the description, serve to explain the principles of the invention: 

[0037] FIG. 1 is a block diagram of a vehicle suspension system according to an 

embodiment of the present invention; 

[0038] FIG. 2 illustrates an exemplary dynamic model for decoupled modal 

equations; 

[0039] FIG. 3 illustrates an exemplary dynamic model for a vehicle suspension 

system including a sky-hook modal damper that satisfies non-coupled modal equations 
according to an embodiment of the present invention; 

[0040] FIG. 4 is a flowchart showing a method for controlling a vehicle 

suspension system, performed by a controller when a number n of degrees of freedom 
of equations of motion of the suspension system equals a number p of actuators, 
according to an embodiment of the present invention; 

[0041] FIG. 5 illustrates an exemplary dynamic model for a vehicle suspension 

system including a sky-hook coulomb modal friction damper that satisfies non-coupled 
modal equations according to an embodiment of the present invention; and 
[0042] FIG. 6 is a flowchart showing a method for controlling a vehicle 

suspension system, performed by a controller when a number n of degrees of freedom 
of equations of motion of the suspension system is greater than a number p of actuators, 
according to an embodiment of the present invention. 
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DETAILED DESCRIPTION OF THE PREFERRED EMBODIMENTS 

[0043] A preferred embodiment of the present invention will hereinafter be 

described in detail with reference to the accompanying drawings. 
[0044] The present invention may be applied to a vehicle regardless of its 

number of wheels (and consequently its number of springs and dampers). However, for 
illustrational purposes, FIG. 1 shows a vehicle suspension system in which four wheels 
120 are included and a corresponding number of springs 140 and dampers 130 are 
mounted to a vehicle body 100. 

[0045] A sensor unit 110 is mounted to the vehicle body 100 so as to measure 

motion parameters of the vehicle body 100. 

[0046] A motion parameter includes each element in a state vector 

x(t) 

(preferably, its time derivative by which an equation of motion of the vehicle 

body is established. The elements of the state vector may be preset in various 

manners according to analysis techniques of rigid body motion. As a preferred example, 

the elements of the state vector 

x(0 

include six parameters for translational and 
rotational motion of the vehicle body 100. The order of the state vector measured by the 
sensor unit 110 relates to degrees of freedom in simultaneous equations of motion of the 
vehicle suspension system. 

[0047] A vehicle suspension system of an embodiment of the present invention 

includes front-left (FL), front-right (FR), rear-left (RL), and rear-right (RR) suspensions. 
[0048] An actuator 160 is mounted to each of the FL, FR, RL, and RR 

suspensions so as to generate and apply activating force to a corresponding suspension. 
FIG. 1 shows four actuators 140, however, the number of the actuators is not required to 
be equal to the number of the wheels 120. 

[0049] A method for designing a vehicle suspension system of a present 

embodiment and a suspension system of a present embodiment are hereinafter described 
in more detail. For the purpose of better understanding the wide scope of possible 
applications of the spirit of the present invention, the number of degrees of freedom in 
an equation of motion of the suspension system is denoted by n, and the number of 
independent actuators is denoted by p. 
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[0050] An equation of motion of vehicle body 100 according to a linear 

suspension system showing viscous damping may be written as the following equation 
(1). 

(1) Mx{t) + C(x(0 - u(tj) + K(x(t) - = Pf{t) 

[0051] Here, M, C, and K respectively denote a mass matrix, a damping 

matrix, and a stiffness matrix, each of which is symmetrically nxn. The mass matrix M 
is a positive definite matrix, the damping matrix C is a positive semi-definite matrix, 
and the stiffness matrix K is a positive definite matrix. 

[0052] The P denotes an nxp real matrix corresponding to positions of the 

actuators, and the denotes a pxl external force (i.e., actuating force generated by 

the actuator 160) vector. 

[0053] *(0and M (0 respectively denote nxl state and disturbance vectors. 

u(t) 

implies disturbance caused while the wheel 120 moves along a road surface. 
[0054] Driving an equation of motion for a given suspension system in the form 

of the above equation 1 is within the capability of a person of ordinary skill in the art. In 

addition, given specific values of the matrices M,C,K , and jP of equation (1), 
extracting physical features of the suspension system corresponding to equation (1) is 
also within the capability of a person of ordinary skill in the art. 
[0055] A decoupling process of equation (1) is firstly described hereinafter. 

[0056] A matrix S that satisfies the following equation (2) can be obtained by 

obtaining eigenvectors of the stiffness matrix K, and subsequently normalizing the 
eigenvector with respect to the mass matrix M. 

(2 ) S T MS = I and S T KS = diag[(D? ] 

[0057] Here, ^iag[cD j ] d enotes a diagonal matrix whose (i,i) element is ®i 

and I denotes an nxn unit matrix. Hereinafter, ^Vk is defined as S T KS (i.e., 
= S T KS y 
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2 

[0058] ®i , being a square of a natural frequency CO, is an i-th eigenvalue of the 

stiffness matrix K. 

[0059] Equation (2) defines an orthonormality of the matrix S. Equation (2) 

implies that, according to a similarity transform by a similarity transform matrix S, the 

mass matrix M may be transformed to a unit matrix, and the stiffness matrix K may 
be transformed to a diagonal matrix. 

[0060] According to the similarity transform by the similarity transform matrix 

S, a similar matrix C (called a modal damping matrix hereinafter) of the damping 
matrix C can be obtained as the following equation (3): 

(3) S T CS = C 

[0061] The modal damping matrix C is an nxn symmetric positive semi- 

definite matrix, and is not necessarily a diagonal matrix. 

[0062] By multiplying the similarity transform matrix S to the left of equation 

(1), and by substitutions of *(0 = ^#00 and «(0 - Sf]{t) equation (1) is 
converted into the following equation (4): 
(4) 

(S T MS)£(t) + (S T CS)(i(t) - 77(0)+ (S T KS)(<*(t) - 77(0) = S T Pf(t) 

[0063] Therefore, if the modal damping matrix can be 

diagonalized, the left side of equation (4) becomes a diagonalized matrix equation, i.e., 
modal equations, by a similarity transform. 

[0064] Such a modal damping matrix C can be diagonalized in various 

manners. As an example, a condition of the following equation (5) may be imposed to 

the suspension system: 

k . = OC XC . 
(5) J J (j=l,..., number of springs) 

SY k C 

[0065] Here, w is a constant, and J and J respectively denote a j-th spring 

coefficient (or equivalently, spring rate) and a j-th damping coefficient. 
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[0066] For a vehicle suspension system, a spring for controlling motion of the 

vehicle body and a damper for damping a vibration of the spring are usually mounted at 

a same position. Therefore, the stiffness matrix K and the damping matrix C show 
same values related to positions, although spring coefficients and damping coefficients 
thereof may be different. 

[0067] Therefore, equation (5) implies that a vehicle suspension system is 

k = a xc . 

designed such that J J is satisfied for every j=l,. . ., number of springs. 

That is, a ratio of spring and damping coefficients of a spring and a damper is the same 
regardless of the position (such as FL, FR, RL, and RR) of the suspension. 

[0068] According to the condition of equation (5), the modal damping matrix C 

A = Sf T KS! 

can be diagonalized, in the same way as the similar matrix K — ° of the 

/\ 

stiffness matrix K . The modal damping matrix C diagonalized as such can be written 
as the following equation (6): 

(6) S T CS = C = diag\2Z i a> i \ 

[0069] Here, £i denotes an i-th modal damping ratio. 

[0070] Rewriting equation 1 (or equation 4) under the condition of equation 5, 

the following equation (7) can be derived: 
(7) 

/|(o + dia&SM ii(t) - m) + a* (m - lit)) = h) 

[0071] Here, / = Qfi* ) and Q = S T P . Regarding X (/ ) and U (0, 

X(0 = -S«0 and "(0 = ^7(0. 

[0072] As is clear from equation (7), the left side of the equation of motion of a 

vehicle suspension system is fully diagonalized. 

[0073] Therefore, n (i=l,. . . ,n) modal equations can be derived from the above 

equation (7) in a form of the following equation (8). 
(8) 
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H0+ 2<r> ( (#, (o - fjj w)+ «>, 2 (?, (o - i7i (0) = / (o = £ Q v fj (o 

[0074] For reference, an exemplary dynamic model for decoupled modal 

equations is illustrated in FIG. 2. 

[0075] As can be seen from equation (8), the equation of motion of a suspension 

system is fully decoupled except for the right side which relates to forced activation. 
And now, the suspension system becomes convenient for designing an actuating force 
of the actuators 160 for vibration control, and/or applying the designed actuating force 
to the suspension system. 

[0076] The above described method for designing a suspension system can be 

summarized as: formalizing the vehicle suspension system by equation (1); calculating 

eigenvectors of a stiffness matrix K of equation (1); normalizing the eigenvectors with 
respect to a mass matrix M of equation (1); calculating a similarity transform matrix o 
consisting of the normalized eigenvectors; and normalizing equation (1) using the 
similarity transform matrix S . 

[0077] According to such a method for designing a vehicle suspension system, 

an equation of motion of a suspension system is mode normalized at the stage of 
designing. Therefore, dynamic behavior characteristics of such a suspension system can 
be easily enhanced by better control logic which can be obtained by an analysis study of 
a vehicle dynamics. 

[0078] Determining active damping of suspension motion embodiments for 

when and how large an actuating force should be generated by the actuators 160 is now 
described in detail, with respect to cases of different numbers of activators. 
[0079] In order for active damping of the suspension motion to occur, a 

suspension system of a present embodiment further includes, as shown in FIG. 1, a 
controller 150 for receiving signals from the sensor unit 1 10 and accordingly for 
controlling the actuators 160 on the basis thereof. 

[0080] The controller 150 can be realized by one or more processors activated 

by predetermined software, and the predetermined software can be programmed to 
perform each step of a method according to a preferred embodiment of this invention. 
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[0081] Hereinafter, a first embodiment is described with respect to a case in 

which the actuators 160 are provided in a same number n of the modal equations (i.e., 
modes), and a second embodiment is described with respect to a case in which the 
number p of the actuators 160 is less than the number n of the modes. 
[0082] Firstly, the former case is described in detail. In this case, the number p 

of the actuators 160 is equal to the number n of modes. In this case, a sky-hook modal 
damper is preferably adopted as a damper for active damping. An exemplary dynamic 
model for a vehicle suspension system including a sky-hook modal damper that satisfies 
non-coupled modal equations is illustrated in FIG. 3. 

[0083] Whereas a conventional viscous damping damper forms damping force 

proportional to relative speed between a vibrating mass (i.e., vehicle body) and a 
supporting plane (i.e., road surface), the sky-hook damper forms a damping force 
proportional to an absolute velocity of the mass (i.e., regardless of undulation of the 
road surface). 

[0084] Therefore, an active damping force (i.e., a modal control force) fi 

designed to such a sky-hook damper may be expressed as the following equation (9): 



(9) 



fi ~ ~C-Si%i 



c 

[0085] Here, Sl is a damping coefficient of a sky-hook damper connected to 

an i-th mode, and %i is an absolute velocity of the i-th mode. 

[0086] When the number p of the actuators 160 equals the number n of the 

modes (n=p), the matrix Q is a regular matrix (or a non-singular matrix) that is 

invertible, and therefore, the actuating force f can be written as the following equation 
(10): 

do) f = Q~ 1 f 

[0087] Using the above equation (10), the equation of motion of the suspension 

system can be converted to the following equation (1 1): 
(11) 

I (O + (2£>, + C Si )g t (0 + a>?& (0 = C.,77, (0 + afrjt (r) 
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[0088] As can be seen from equation (1 1), the equation of motion of a 

suspension system is fully decoupled. Therefore, when an actuating force of equation 
(10) is applied to the actuators 160, a modal vibration in each mode of the equation of 
motion of the suspension system can be easily diminished. 

[0089] The controller 150 can apply the actuating force of equation (10) to the 

actuators 160 in the process described hereinafter. 

[0090] As shown in FIG. 4, firstly at step S410, the controller 150 calculates the 

time derivative (i.e., velocity) of the state vector on the basis of signals 

from the sensor unit 110. 

[0091] Subsequently at step S420, the controller 150 calculates a velocity 

of a modal state vector on the basis of the velocity of the state vector 

x(r) 

The velocity #(0 of the modal state vector #(0 is calculated by 

using the similarity transform matrix S and the stiffness matrix K , which is 

predetermined at a stage of designing the suspension system. In more detail, the velocity 

of the modal state vector is calculated as a value according to 

i(t) = (S T KS)-' (S T K)x(t) 

[0092] Subsequently at step S430, the controller 150 calculates an actuating 

force vector for the actuators 160 using the velocity of the modal state 

vector 

[0093] In more detail, the actuating force vector is calculated by an 

equation 

f(t) = {s , P)-\-c Si )4(t) 

equivalent to equation (10) using the 
similarity transform matrix S and the matrix P corresponding to positions of the 
actuators 160, which is predetermined at a stage of designing the suspension system. 

The actuating force vector calculated as such is calculated as a value according 

t0 m = (,s T py (-c sl )(s T Ksy (s r *)*(o. 

[0094] Subsequently at step S440, the controller 150 applies such calculated 

actuating force vector to the actuators 160. 
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[0095] Regarding the above description, the following remarks are notable. 

[0096] Firstly, it is not necessarily required that the sensor unit 1 10 detects the 

state vector and the controller 150 calculates its velocity . As a modified 

embodiment, the sensor unit 1 10 may directly detect the velocity of the state 

vector such that the controller 150 can simply receive the detected velocity 

Therefore, in the detailed description of the embodiment of the present invention and in 

the appended claims, the calculation of the velocity of the state vector 

x(t) 

should be understood to cover any step for obtaining the value of the velocity of 
the state vector 

[0097] In addition, a process for calculating the actuating force vector 

AO 

has been described with reference to a plurality of steps in the above description, 
however it should be understood to be an exemplary process to finally calculate a value 

of f{t) = (S T Py (-C Si )(S T KSr (S T K)m. As a s, m ple variation, 

having stored a value of in the controller 150, 

the controller 150 may obtain the actuating force vector by simply multiplying 

the stored matrix by the velocity vector 

after the step S410. 

[0098] A suspension system according to a second embodiment of the present 

invention is hereinafter described in detail. 

[0099] According to a second embodiment, the number p of the actuators 160 is 

less than the number n of modes (or equivalently, the number of modal equations). 
[00100] In this case, a sky-hook coulomb modal friction damper is preferably 

adopted as a damper for active damping. An exemplary dynamic model for a vehicle 
suspension system including a sky-hook coulomb modal friction damper that satisfies 
non-coupled modal equations is illustrated in FIG. 5. 

[00101] When the number p of the actuators 160 is less than the number n of the 

modes, the equation of motion of the suspension system can not be entirely decoupled 
including its right side (i.e., actuating force side). 
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[00102] In this case, when an actuating force proportional to a mode velocity is 

applied to the actuators 160 in order to increase a value of a diagonal element of the 

modal damping matrix C , values of its non-diagonal elements are also increased. That 
is, damping effect cannot be independently increased with respect to a specific mode, 
and an increase of damping for a specific mode may cause a negative effect (i.e., 
increase of vibration energy) on damping in a different mode. 

[00103] Therefore in such a case in which the number p of the actuators 160 is 

less than the number n of the modes, a sky-hook coulomb modal friction damper is 
adopted as a damper for an active damping such that a negative effect on different 
modes when damping of a specific mode is controlled to increase. 
[00104] In the case of a sky-hook coulomb modal friction damper, the modal 

control force f satisfies the following equation (12). 
(12) 

/, = -F Si sign(l) = £&,./, 

[00105] Here, Si denotes a frictional force of a sky-hook coulomb friction 

damper connected to an i-th mode. As was previously described, p denotes the number 

of the actuators, and the matrix Q is defined as Q ~~ ^ ^* . means an absolute 

velocity of an i-th mode, and signed denotes a sign thereof. 

f. 

[00106] In order to satisfy the above equation (12), a j-th actuating force J J 

(j=l,...,p) is calculated as the following equation (13): 
(13) 

V0^igw(£)^0&a//gn(fe)>0&..^/gn(^)^Q, fj=-F p 



\ifQ/ig<&<0&Q v sign&)<0&. • Qfig^O <<X fj =-F N 
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F F 

[00107] Here, P and N respectively denote a positive maximum force and a 



negative maximum force that a j-th actuator can generate 

F F 
is a value between P and N m 



F k(ot k = \,...{2"-2) 



F 

[00108] According to the above equation (13), when r k is not zero (0) for 



k=l,...{2 n -2) 



motion of the vehicle body 100 is reduced for several modes and 

increased for several other modes. For example, in order to reduce motion of the first 
i? 

mode, 1 should have a positive value, but in this case, the motion increases in an n-th 
mode. Therefore, according to a second embodiment, values k for all 

Jr — 1 (2* —2) 

9 ' ' ' >\ / are set to zero (0) such that motions in each mode can be 

independently reduced. 

f 

[00109] In summary, the controller 150 calculates the j-th actuating force J J 

(j=l,...,p) as the following equation (14): 
(14) 

V^^(§)^0>&ri=l,...A fj=~F P 
ifQySign&KO jbri=\.,./i> fj=-F N 
Otherwise 9 f . = 0 

V J 

[00110] According to the above equation (14), the actuators 160 are activated 

only in the case in which vibration energy can be dissipated in every mode. 
[00111] For example, when only two actuators are mounted to a vehicle 

suspension system that has more than two degrees of freedom, only a first actuator is 
activated with an actuating force if the first actuator dissipates vibration energy in every 
mode but a second actuator does not. In other words, a j-th actuator is operated only 
when the operation ensures dissipation of vibration energy in every mode. 
[00112] A first line of equation (14) (or a first line of equation (13)) implies that, 

in the case in which the j-th actuator can dissipate vibration energy in every mode, the j- 
th actuator is operated to the positive maximum force. That is, in this case application of 
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a positive actuating force to the j-th actuator dampens vibration. So, such vibration 
damping effect is controlled to be maximized by applying the positive maximum 
actuating force to the j-th actuator. 

[00113] In addition, a second line of equation (14) (or a last line of equation (13)) 

implies a case in which a positive value of an actuating force of the j-th actuator causes 
an increase of vibration energy in every mode. So, in this case, the j-th actuator is 
controlled to a negative maximum actuating force. That is, since application of a 
positive actuating force to the j-th actuator increases vibration, such vibration can be 
dampened by applying negative actuating force. So, in this case, vibration damping 
effect is controlled to be maximized by applying the negative maximum actuating force 
to the j-th actuator. 

[00114] Lines between the first and last lines of equation (13) relate to cases in 

which operation of the j-th actuator causes an increase of vibration for several modes 
and a decrease of vibration for several other modes. A person of ordinary skill in the art 

may choose appropriate values of r k to comply with a desired characteristic of a 
specific suspension system within the prescribed ranges. 

[00115] According to the second embodiment, the values of 17 k are preset to zero 

(0) such that any negative effect (increase of vibration energy in any mode) may not 
occur. 

[00116] The controller 150 can apply the actuating force of equation (14) to the 

actuators 160 in the process described hereinafter. 

[00117] As shown in FIG. 6, firstly at step S610, the controller 150 calculates the 

time derivative (i.e., velocity) of the state vector on the basis of signals 

from the sensor unit 1 10. 

[00118] Subsequently at step S620, the controller 150 calculates a velocity <f(0 

of a modal state vector on the basis of the velocity of the state vector 

The velocity 4(f) of the modal state vector #(0 is calculated by 

using the similarity transform matrix S and the stiffness matrix K , which is 

predetermined at a stage of designing the suspension system. In more detail, the velocity 
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<f(0 of the modal state vector is calculated as a value according to 

i(t) = (S T KSy l (S T K)x(t\ 

[00119] Subsequently at steps S630-S640, the controller 150 calculates the j-th 

f. 

component J J (i.e., an actuating force of a j-th actuator) of the actuating force vector 
f(t) for each j (j=l,...,p), using the velocity of the modal state vector 

4(f). 

f. 

[00120] In order to calculate the actuating force J J for the j-th actuator, the 

controller 150 firstly determines if Qtj^S n (^i) — 0 j s satisfied for every i-th mode 
(i=l,...,n) at step S630. 

[00121] When Qy Si S n (^i) ^ 0 is satisfied for every i-th mode (i=l,...,n) 

f -F 
(S630-yes), the controller 150 sets the actuating force J J for the j-th actuator as 1 P at 

step S635. 

[00122] In addition, the controller 150 further determines if Qv si S n (^i) < 0 j s 

satisfied for every i-th mode (i=l,...,n) at step S640. 

[00123] When Qv si S n i^i) < 0 is satisfied for every i-th mode (i=l,...,n) 

f —F 

(S640-yes), the controller 150 sets the actuating force J J for the j-th actuator as N 
at step S645. 

[00124] When neither of Qv si S n i^i) ^ 0 and QijSlgni^d < 0 ig fuUy 

satisfied for all i-th modes (i=l,...,n) (S630-no and S640-no), the controller 150 sets the 
f 

actuating force J J zero (0) at step S650. 

[00125] By recursive execution of the steps S630-S650 according to step S660, 

f 

the controller 150 calculates the actuating force J J for the j-th actuator (j=l,...,p). 

f. 

[00126] When the actuating force J J has been calculated for every j=l,. . .,p, the 

calculation of the actuating force vector AO is completed. 

[00127] Then at step S670, the controller 150 activates the actuators 160 by the 

calculated actuating force vector 
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[00128] Regarding the above description, the following remarks are noteworthy. 

[00129] Firstly, it is not necessarily required that the sensor unit 1 10 detects the 

state vector and that the controller 150 calculates its velocity .Asa 

modified embodiment, the sensor unit 1 10 may directly detect the velocity of the 

state vector such that the controller 150 can simply receive the detected velocity 

. Therefore, in the detailed description of the embodiment of the present invention 
and in the appended claims, the calculation of the velocity of the state vector 

should be understood to cover any step for obtaining the value of the velocity 
m of the state vector 

[00130] According to an embodiment of the present invention, an actively 

controlled suspension system can be more easily designed for its active control. 
[00131] When as many actuators are provided as there are degrees of freedom of 

the suspension system, each vibration mode can be efficiently and independently 
controlled. 

[00132] In addition, vibration modes can be controlled without causing negative 

effects to other vibration modes even when the number of actuators falls short of the 
number of degrees of freedom of the suspension system. 

[00133] While this invention has been described in connection with what is 

presently considered to be the most practical and preferred embodiment, it is to be 
understood that the invention is not limited to the disclosed embodiments, but, on the 
contrary, is intended to cover various modifications and equivalent arrangements 
included within the spirit and scope of the appended claims. 
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